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Abstract. We present an updated phase-coherent timing solution for 
the young, energetic pulsar PSR B1509— 58 for twenty years of data. 
Using a partially phase-coherent timing analysis, we show that the sec- 
ond frequency derivative is changing in time implying a third frequency 
derivative of v= (—9 db 1) X lCP 32 s -4 . This value is consistent with the 
simple power law model of pulsar rotation. 



1. Introduction 

Radio pulsars are powered by rotational kinetic energy and experience losses due 
to the emission of electromagnetic energy while spinning down. The spin-down 
of radio pulsars is described by (Manchester & Taylor, 1977) 

z>oc-i/\ (1) 

where v is the spin frequency, i> is the frequency derivative and n, the braking 
index, equals 3 for a spinning magnetic dipole in a vacuum. Taking a derivative 
of (1), we find that n = vvjv 2 . Taking further derivatives of (1), and defining 
the second braking index as tuq = n(2n — 1), we find an expression m = v 2 v jifi . 
Measuring a value of m that agrees with the predicted value mo is an excellent 
check on the validity of the power law spin-down model of radio pulsars given 
in equation (1). Kaspi et al. (1994) reported a value of n = 2.837 ± 0.001 and 
m = 14.5 ± 3.6 for PSR B1509-58. 



2. Results and Discussion 

Performing a phase-coherent timing analysis on twenty years of radio timing 
data for PSR B1509— 58 obtained from the Molongolo Observatory Synthesis 
Telescope and the Parkes 64m telescope, we found the second frequency deriva- 
tive to be v = 1.962(1) x 10~ 21 s -3 , implying a braking index of n = 2.814(1). 
Significant low-frequency timing noise in the data prohibited a phase-coherent 
analysis to determine an accurate value of the third frequency derivative. This 
prompted a partially phase-coherent timing analysis by breaking up the data 
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Figure 1. Phase-coherent subsets plotted versus epoch. The slope of 
the line is the third derivative and has a value of (—9 ± 1) x 10~ 32 s -4 . 



into subsets of ~ 2 yr (such that phase residuals were 'white' in each interval 
following a fit of v, i> and v) and plotting the values of v against time, the results 
of which are shown in Figure 1. This method is less sensitive to timing noise 
which dominates on long time scales. We found a value v= (— 9± 1) x 10 _32 s 
by a weighted least-squares fit to the data. This value for v is in agreement with 
that predicted from the spin-down law of z>= —9.37 x 10~ 32 s~ 4 . This measured 
value for i) implies a second braking index of m = 12.7 ± 1.4, also in agreement 
with the predicted value of m = 13.3, though less than the value m = 15 ex- 
pected for an ideal spinning dipole. Our measurement of the v implies that the 
spin-down law correctly predicts the behaviour of PSR B1509— 58, reinforcing 
the use of the power law given in equation (1) to describe the rotation of pulsars. 
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